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Abstract 

We develop a theory of existence and uniqueness for the following porous 
medium equation with fractional diffusion, 

^ + (-Ay/ 2 (\u\ m ~ 1 u) = 0, xeR N ,t>0, 
at 

u(x,0) = f(x), xeR N . 

We consider data / G L 1 (M Ar ) and all exponents < a < 2 and m > 0. 
Existence and uniqueness of a weak solution is established for m > m* = (N — 
cr) + /N, giving rise to an /^-contraction semigroup. In addition, we obtain the 
main qualitative properties of these solutions. In the lower range < m < m* 
existence and uniqueness of solutions with good properties happen under some 
restrictions, and the properties are different from the case above m*. We also 
study the dependence of solutions on /, m and a. Moreover, we consider the 
above questions for the problem posed in a bounded domain. 
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1 Introduction 



The aim of this paper is to develop a theory of existence and uniqueness, as well as to 
obtain the main qualitative properties, for a family of nonlinear fractional diffusion 
equations of porous medium type. More specifically, we consider the Cauchy problem 

— + (-A) a / 2 (\u\ m - l u) = 0, x£R N ,t>0, 
dt 

u(x,0) = f(x), xeR N . 

We take initial data / G L 1 ^), which is a standard assumption in diffusion prob- 
lems, with no sign restriction. As for the exponents, we consider the fractional expo- 
nent range < a < 2, and take porous medium exponent m > 0. In the limit a — > 2 
we recover the standard Porous Medium Equation (PME) 

which is a basic model for nonlinear and degenerate diffusion, having now a well- 
established theory [39]. 

The nonlocal operator (—A) "/ 2 , known as the Laplacian of order a, is defined for 
any function g in the Schwartz class through the Fourier transform: if (— A) a / 2 g = h, 
then 

(i-2) Mo = ie?(o- 

If < a < 2, we can also use the representation by means of an hypersingular kernel, 

(1-3) (-A)^) = C N , a P.V. / 9{X) ~ |g dz, 

J r n \x-z\ N+a 

where Cn,o = 2 ^n/^i^J^ * s a norma li za tion constant, see for example [30]. There 
is another classical way of defining the fractional powers of a linear self-adjoint non- 
negative operator, in terms of the associated semigroup, which in our case reads 

(1-4) ("Ar /2 s(*) = --^ \ {e tA g(x) - g(x)) ™. 

It is easy to check that the symbol of this operator is again The advantage of 
this approach is that it gives a natural way of defining the problem in a bounded 
domain, by means of the spectral characterization of the semigroup e tA , see [SB]- The 
problem posed in a bounded domain is also studied in this paper. 

Equations of the form (II. ip can be seen as fractional- diffusion versions of the PME. 
Though our paper is aimed at providing a sound mathematical theory for this evo- 
lution equation and the nonlinear semigroups generated by Problem (11. ip and the 
problem posed on bounded domains, we mention that the equation appears as a 
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model in statistical mechanics and the linear counterpart in [27]. We also want 
to point out that there are other natural options of nonlinear, possibly degenerate 
fractional-diffusion evolutions under current investigation. Thus, the papers [H], [T5] 
consider the following fractional diffusion PME dtu — V • (uV(— A)~ s ^ 2 u) = 0. It has 
very different properties from the ones we derive for (11.11) . The standard PME (with 
m = 2) is recovered in such model for s = 0. A more detailed discussion on these 
issues is contained in the survey paper [4TJ] . 

These two kinds of equations can be also viewed as nonlinear versions of the linear 
fractional diffusion equation obtained for m = 1, which has the integral representation 



where K a has Fourier transform K a (^,t) = e~'^ CT *. This means that, for < a < 2, 
the kernel K a has the form K a (x,t) = t~ N ' a F(\x\ t~ l l a ) for some profile function F 
that is positive and decreasing, and behaves at infinity like F{r) ~ r -{N+a) jgj When 
a = 1, F is explicit; if a = 2 the function K 2 is the Gaussian heat kernel. The linear 
model has been well studied by probabilists, since the fractional Laplacians of order 
a are infinitesimal generators of stable Levy processes [I], [7J. However, an integral 
representation of the evolution like (jl.5p is not available in the nonlinear case, thus 
motivating our work. 

In a previous article [35] we studied Problem (II. ip for the particular case a = 1. The 
key tool there was the well-known representation of the half-Laplacian in terms of 
the Dirichlet-Neumann operator, which allowed us to formulate the nonlocal problem 
in terms of a local one (i.e., involving only derivatives and not integral operators). 
For a ^ 1, Caffarelli and Silvestre p3] have recently given a similar characterization 
of the Laplacian of order a in terms of the so-called a-harmonic extension, which is 
the solution of an elliptic problem with a degenerate or singular weight. However, 
even with this local characterization at hand, many of the proofs that we gave for 
cr = 1 cannot be adapted to cope with a general a. Hence we have needed to use new 
tools, which in several cases do not involve the extension. These techniques, which 
include some new functional inequalities, have also allowed us to improve the results 
obtained in [35] for the case a = 1. The case of a bounded domain was not treated 
in that paper. 

2 Main results 

As in the case of the PME, there is a unified theory of existence and uniqueness 
of a suitable concept of weak solution above a critical exponent, given for a general 
a G (0, 2) by m, = (N — a) + /N. The linear case m = 1 lies always in this range. 

Theorem 2.1 Let m > m* and a G (0,2). For every / G L 1 (IR iV ) there exists a 
unique weak solution of Problem ( II. ip . 
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Figure 1: The critical line m* = (N — o~)/N and the supercritical region m > m* for 
N > 3, N = 2, and N = 1. 

The precise definition of weak solution that guarantees uniqueness is stated in Def- 
inition 13.11 or equivalently in Definition 13.21 

The construction of the solution in the previous theorem follows from a double limit 
procedure, approximating first the initial datum by a sequence of bounded functions, 
and also approximating M. N by a sequence of bounded domains with null boundary 
data. In this respect we show existence of a weak solution to the associated Cauchy- 
Dirichlet problem, a result that has an independent interest, see Section HI 

The weak solutions to Problem f ll.ip have some nice qualitative properties that are 
summarized as follows. 

Theorem 2.2 Assume the hypotheses of Theorem \2.\\ and let u be the weak solution 
to Problem f ll.ip . 

(i) d t u G L°°((r, oo) : L 1 (R N )) for every r > 0. 

(ii) Mass is conserved: J RN u(x, t) dx = f RN f(x) dx for all t > 0. 

(iii) Let Mi,«2 be the weak solutions to Problem fll.ip with initial data /i,/2 G 
L 1 (R JV ). Then, 

/ {u\ — u 2 ) + {x,i) dx < l (fi — f2)+{x) dx, (L 1 -order-contraction property). 
Jr n Jr n 

(iv) Any LP -norm of the solution, 1 < p < oo, is non-increasing in time. 

(v) The solution is bounded in WL N x [r, oo) for every r > ( L 1 -L°° smoothing 
effect/ Moreover, for all p > 1, 

(2-1) ll«(-,*)llL-(a*)< ^"^11/11^). 

with 7 P = (m — 1 + ap/N)~ l , 5 P = ap'jp/N, and C = C(m,p, N, a). 

(vi) If f > the solution is positive for all x and all positive times. 
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(vii) If either m > 1 or f > 0, then u G C a (M. N x (0, oo)) for some < a < 1. 

(viii) TTie solution depends continuously on the parameters a G (0,2), m > m* ; and 
f G L 1 (R JV ) in the norm of the space C([0,oo) : L X (R^)). 

Remarks. (a) Properties (i) and (ii) were only known for a = 1 in the case of 
nonnegative initial data [35] . 

(b) The positivity property (vi) is not true for the PME in the range m > 1. The 
fact that it holds for Problem (II. ID stems from the nonlocal character of the diffusion 
operator. 

(c) A weak solution satisfying property (i) is said to be a strong solution, see Defini- 
tion 16.11 These kind of solutions satisfy the equation in (11.11) almost everywhere in 
Q = R N x (0,oo). 

Our main interest in this paper is in describing the theory in the above mentioned 
range m > m*. However, we also consider the lower range < m < m* for contrast. 
In that range (which implies that 0<a<lifA^=l)we obtain existence if we 
restrict the data class (or if we relax the concept of solution). In addition, in order 
to have uniqueness, we need to ask the solution to be strong. 

Theorem 2.3 Let a G (0,2) and < m < m*. For every f G L 1 (R JV ) n L P (R N ) 
with p > p*(m) = (1 — m)N/a there exists a unique strong solution to Problem (11.11) . 




Figure 2: The critical line = (1 — m)N/a and the existence region of strong 
solutions. 

This theorem improves the results obtained of [53] for a — 1, since in that paper 
existence and uniqueness of a (weak) solution were only proved for the case of inte- 
grable and bounded initial data, / G L X (R N ) nL°°(M Ar ). Moreover, the weak solution 
was only proved to be strong in the case of nonnegative initial data. 

Some of the properties of the solutions in this lower range are rather different from 
those in the upper range. 
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Theorem 2.4 Assume the hypotheses of Theorem \'2.'d\ and let u be the strong solu- 
tion to Problem (11.11) . 

(i) The mass J RN u(x,t)dx is conserved if m = m*. Mass is not conserved if 
m < m* . Actually, when < m < m* there is a finite time T > such that 
u(x,T) = inR N . 

(ii) There is an L 1 - order- contraction property. 

(iii) Any LP -norm of the solution, 1 < p < oo, is non-increasing in time. 

(iv) The solution is bounded in R N x [r, oo) for every r > (L p -L°° smoothing 
effect). Moreover, if p > p*(m) (which is necessary to make ^ p > 0) then 
formula ( 12.11) holds. 

(v) If f > the solution is positive for all x and all positive times up to the extinc- 
tion time. 

(vi) Let f > and let T be the extinction time. Then u G C a (M. N x (0, T)) for some 
< a < 1. 

We will make a series of comments on these two results. First, we point out that 
the result on the conservation of mass for m = m* is new even for a = 1. 

A more essential observation is that there is an alternative approach: using the 
results of Crandall and Pierre [22] we can obtain the existence of a unique so-called 
mild solution for all / G L 1 (IR Ar ) in the whole range of m and a, via the abstract 
theory of accretive operators. This approach has therefore the advantage of having 
general scope. Two problems arise with this way of looking at the problem: (a) how 
to characterize the mild solution in differential terms; and (b) how to derive its 
properties. Our paper answers both questions. 

The strong solutions that we have constructed are mild solutions. Hence, in our 
restricted range of initial data the unique mild solution is a strong solution. For a 
general / G L 1 (M iV ) and m < m„ we will show that the mild solution is a very weak 
solution (a solution in the sense of distributions), see Theorem 18.41 However, we will 
fail to prove that this very weak solution is a weak solution in the sense of Definition 
13. 1[ and hence we will not be able to obtain the properties listed above (bit note that 
the L 1 -contraction holds since it is a consequence of the accretivity of the operator). 

As to the continuous dependence of the solution in terms of the parameters, conver- 
gence in L 1 (M Ar ) is not expected to hold for < m < m*, since mass is not conserved 
in that region. Instead we expect to have continuity in weighted spaces, much in the 
spirit of [5]. Nevertheless, we are able to extend the continuous dependence result of 
Theorem 12.21 (viii) to cover the case m = for N > 2, see Proposition llO.il We also 
show that continuity holds in the upper limit a — > 2, thus recovering the standard 
PME, see Theorem fT02l 
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3 Weak solutions. An equivalent problem 



3.1 Weak solutions 

If ip and tp belong to the Schwartz class, Definition (II. 2p of the fractional Laplacian 
together with Plancherel's theorem yield 
(3.1)^ 

(-A)*/ 2 ^= / \£\*$<p= I \^ /2 m a/2 0= I (-A)*/V (-A) ff /V 



Therefore, if we multiply the equation in (II -ip by a test function p> and integrate by 
parts, as usual, we obtain 

(3.2) r f u^fdxds - [°° f {-^Y l \\u\ m ~ 1 u){-^Y l ^dxds = Q. 

Jo Jr n Jo Jr n 

This identity will be the basis of our definition of a weak solution. 

The integrals in (I3.2p make sense if u and |w| m_1, u belong to suitable spaces. The 
correct space for |M| m_1 M is the fractional Sobolev space H a / 2 (R N ), defined as the 
completion of C^(R N ) with the norm 

imIjw>= (7 ieri^r^y /2 = ii(-Ar/vii 2 . 

\Jr n / 

Definition 3.1 A function u is a weak (L l -energy solution) to Problem (II. ip if: 

• u G C([0,oo) : L\R N )), \u\ m ^u G ^((0,00) : H^ 2 (R N )); 

• identity (ET2|) holds for every tp G Cq(R n x (0, 00)); 

• u(-,0) = f almost everywhere. 

For brevity we will call weak solutions the solutions obtained below according to this 
definition, but the complete name weak ^-energy solution is used in the statement 
to recall that we are making a very definite choice. 

The main disadvantage in using this definition is that there is no formula for the 
fractional Laplacian of a product, or of a composition of functions. Moreover, there is 
no benefit in using compactly supported test functions since their fractional Laplacian 
loses this property. To overcome these difficulties, we will use the fact that our solution 
u is the trace of the solution of a local problem obtained by extending |M| m_1 M to a 
half-space whose boundary is our original space. 

3.2 Extension Method 

If g = g{x) is a smooth bounded function defined in R N , its cr-harmonic extension to 
the upper half-space, v = E(g), is the unique smooth bounded solution v = v(x, y) to 

V ■ (y^Vv) = 0, in R* +1 = {(x, y) G R N+1 : x G R N , y > 0}, 

v(x,0)=g(x), xeR N . 
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(3.3) 



Then, Caffarelli and Silvestre [13] proved that 

(3.4) - ^ lim y 1 -^ = (-Ay' 2 g(x), ^ = 2- 1 I>/2)/r(l - a/2). 

y ^o+ dy 

In (13.31) the operator V acts in all (x, y) variables, while in f !3.4p (—A) "/ 2 acts only 
on the x = (x±, ■ ■ ■ ,Xn) variables. In the sequel we denote 

L a v = V • (y'-'Vv), %L = Pa lim y 1 -^. 

dy a y->o+ dy 

Operators like L a , with a coefficient which belongs to the Muckenhoupt space 
of weights A2 if < a < 2, have been studied by Fabes et al. in [21]. We make use of 
this theory later in the proof of positivity, see Theorem 19.11 

With the above in mind, we rewrite Problem (jl.ip as a quasi- stationary problem for 
w = E(|M| m_1 M) with a dynamical boundary condition 

L a w = 0, (x,y) G R^ +1 , t > 0, 

{ — - l — =0, X6t^ = 0,t>0, 
w = l/l™" 1 /, x G R N , y = 0, t = 0. 

To define a weak solution of this problem we multiply formally the equation in (13.51) 
by a test function ip and integrate by parts to obtain 

f°° [dp f°° f 

(3.6) / / u— dxds — fi a I / y 1 ~ a (V 'w,V '<p) dxdyds = 0, 

Jo Jr n "t J J r n+i 



where it = | Tr(u>)|™ Tr(u>). This holds on the condition that p vanishes for t = 
completion of C^°(M+ +1 ) with the norm 



and also for large y and t. We then introduce the energy space X CT (IR+ +1 ), the 



1 1/2 

(3.7) IMU ct = ( fJ'cr I y 1 ~ a \\7v\ 2 dxdy 



JV+l 



The trace operator is well defined in this space, see below. 
Definition 3.2 A pair of functions (u,w) is a weak solution to Problem (13.5p if: 

• u = I TrHI^ 1 Tr(w) G C([0, 00) : L 1 (M Ar )) J u> G Lg, c ((0, 00) : X CT (M^ +1 )); 

• ideniity (JUS) holds for every ip G Cq 1 x (0, 00) 

• u(-,0) = f almost everywhere. 
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For brevity we will refer sometimes to the solution as only u, or even only w, when no 
confusion arises, since it is clear how to complete the pair from one of the components, 
u = | Tr(w)!^ 1 Tr(V), w = EiH^u). 

The extension operator is well defined in H a / 2 (M. N ). It has an explicit expression 
using a cr-Poisson kernel, and E : H a / 2 (M. N ) — > X a (R+ +1 ) is an isometry, see [13]. The 
trace operator, Tr : X°"(M+ +1 ) — > H a / 2 (R N ) is surjective and continuous. Actually, 
we have the trace embedding 

(3.8) \\<S>\\ X « > || E(Tr($))|| x . = || Tr($)||^ CT/2 

for any $ G X°{R N ). 

3.3 Equivalence of the weak formulations 

The key point of the above discussion is that the definitions of weak solution for our 
original nonlocal problem and for the extended local problem are equivalent. Thus, 
in the sequel we will switch from one formulation to the other whenever this may 
offer some advantage. 

Proposition 3.1 A function u is a weak solution to Problem (II. ip if and only if 
(u, E(|M| m_1 ^)) is a weak solution to Problem (13. 5p . 

Since E : H a/2 (R N ) ->■ X CT (M^ +1 ) is an isometry, we have 

f ^(VEW,VE( V ))=/ (-A) CT /V (-A) CT /V, 

for every ip, ip G H a / 2 (M. N ). Hence the result follows immediately from the next 
lemma, which states that any u-harmonic function is orthogonal in X°"(K+ +1 ) to 
every function with trace on R . 

Lemma 3.1 Letip G H a/2 (R N ) <md$i,$ 2 e -X" CT (R+ +1 ) such that Tr($i) = Tr($ 2 ). 
Then 

/v I y l ~ a (V E(if>), V$i> = \t a [ y'-^V EtyO, V$ 2 ). 

Proof. Let h = $i — $2- Since E(^/>) is smooth for y > 0, given e > we have, after 
integrating by parts, 

y, a ( f y 1 ~ a (VE(il>),Vh)dxdy = ii <r [ £ 1 ~ <T ^^(x, e)h(x, e) dx. 

Je Jr n Jr n oy 

The left-hand side converges to ji a J r n+i y 1_cr (V E(ip), V/i), while the right hand side 
tends to 0, since identity (13 .4p holds in the weak sense in H~ a / 2 (R N ), and Tr(h) = 0. □ 
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4 The problem in a bounded domain 



As an intermediate step in the development of the theory for Problem Q1.1J) . we will 
also consider the Cauchy-Dirichlet problem associated to the fractional PME, 



(4.1) 



7? + (-^) a/2 (\u\ m " l u) = 0, iefi,t>0, 
at 

u = 0, x G <9ft, t > 0, 

k u(x,0) = /(x), x G ft, 



where ft C ffi^ is a smooth bounded domain, / G L 1 (ft). This problem has an interest 
in itself. 

Let us present here the main facts and results about this problem. In view of 
formula (jl.4p . the fractional operator (—A) "/ 2 in a bounded domain can be described 
in terms of a spectral decomposition. Let {ifk, \k}kLi denote an orthonormal basis of 
L 2 (ft) consisting of eigenfunctions of —A in ft with homogeneous Dirichlet boundary 
conditions and their corresponding eigenvalues. The operator (—A) "/ 2 is defined for 
any u G Cg°(fi), u = J2T=i u k l Pk, by 

oo 

(-A)-/ 2 n = ^Af um- 

k=l 

This operator can be extended by density for u in the Hilbert space 

oo 

Hl /2 (Sl) = {ue L 2 (ft) : \\u\\^ /2 = J2K /2 4 < oo}- 

k=l 

Definition 4.1 A function u is a weak solution to Problem ( 14.1 \ if: 
• u G C([0,oo) : L\tt)), lul^u G L? oc ((0,oo) : H° /2 {n)); 



Identity 



[ [ u^fdxds- [ [ (-A) a/4 u m (-A) a/4 ipdxds = 
Jo Jq fit Jo Jq 



holds for every ip G Cq(Q x (0, T)); 
• m(-,0) = / almost everywhere in ft. 

The hypotheses that we need in order to have existence when the spatial domain is 
bounded coincide with the ones we have when the spatial domain is the whole Mr. 

Theorem 4.1 Problem ( 14. ip has a unique weak solution if m > m* and f G L 1 (ft) ; 
which is moreover strong, and a unique strong solution if m < and f G L p (ft) 
with p > p*(m) = (1 — m)N/a. 
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As for the properties of the solutions, most of them, though not all, coincide with the 
ones that hold when the domain is the whole space. 



Theorem 4.2 Assume the hypotheses of Theorem J^.l, and let u be the strong solu- 
tion to Problem (14.11) . 

(i) The solution is bounded in Q x [r, oo) for every r > 0. Moreover, a formula 
analogous to ( 12. ip holds. 

(ii) As a consequence, j n u(x,t) dx = 0(t~ lp ). Moreover, if < m < 1 there is 
extinction in finite time. 

(hi) There is an L 1 -order- contraction property. 

(iv) Any LP -norm of the solution, 1 < p < oo, is non-increasing in time. 

(v) The solution depends continuously on the parameters a G (0,2), m > m*, and 
f G L x (f2) in the norm of the space C([0, oo) : L l (Q)). 

The results of [3] imply that u G C a for m > 1. Positivity for any m > when the 
initial data are nonnegative, and C a regularity for m < 1 are still open problems. 

The construction of a solution uses the analogous to the Caffarelli-Silvestre exten- 
sion (I3.3p . restricted here to the upper half-cylinder Cn = Q x (0, oo), with null 
condition on the lateral boundary, dQ x (0,oo), a construction considered in [T2] . 
[38], pm, [15]. Thus, w = E(|«r -1 u) satisfies 



(4.2) 



Law = 0, (x, y) G Cn, t > 0, 

w = 0, x G <9fi,y > 0, t > 0, 

dy- dt = °> xGfi,y = 0,t>0, 

I ^=|/| m -7, xGfi,y = 0, t = 0. 



In order to define a weak solution to (14. 2 p we need to consider the space Xq (Cn), 
the closure of C£°(Cn) with respect to the norm ( ET7|) . with substituted by C u . 

The extension and trace operators between H^ 2 (Q) and Xq (Cn) satisfy the same 
properties as in the case of the whole space. In fact 



1/2 

'Miff*/* 

1/2 



(*,J c I/'-lVEMI 3 ) =||E(^)|k- 

= ||(-A)-/V||a=(E^iA^ 
See for instance [10] for the explicit expression of E(y?) in terms of the coefficients ipk. 
Definition 4.2 A pair of functions (u,w) is a weak solution to Problem (14. 2 p i/; 
• u=\ Tr(w)\^- 1 Tr(w) G C([0, oo) : L 1 ^)), w G Lf oc ((0, oo) : X CT (Cn))/ 
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• identity 

f°° f dip f°° f 

/ / u— dxds — fi a / / y l ~ a (Vu>,Vy2) dxdyds = 0, 
Jo Jn &t Jo Jc n 

holds for every ip = <p(x, y, t) such that ip G C 1 (Q x [0, oo) x (0, oo)); 

• w(-,0) = / almost everywhere in Q. 

As it happens for the case where Q = R , the two definitions of weak solution, 
Definitions 14.11 and 14.21 are equivalent. 

Remark. The space H^ 2 (Q) can also be defined by interpolation, see [32]. We notice 
that, though for 1 < a < 2 the solutions are zero almost everywhere at the boundary, 
for < a < 1 the functions in H^ 2 {Vt) do not have a trace, [32]- Therefore, the 
boundary condition must be understood in a weak sense, see also [12] for the case 
a = 1. 



5 Some functional inequalities 

In this section we gather some functional inequalities related with the fractional Lapla- 
cian, both in the whole space or in a bounded domain, that will play an important 
role in the sequel. The first one, Strook-Varopoulos' inequality, is well known in the 
theory of sub-Markovian operators [33]. Nevertheless, we give a very short proof using 
the extension operator that makes apparent the power of this technique. 

Lemma 5.1 (Strook-Varopoulos' inequality) Let < 7 < 2 ; q > 1. Then 

(5.1) [ (M'-^X-A)^ > ^Z?± [ |(-A)^ 4 M*/ 2 | 2 

Jrn q J r n 

for all v e L q (R N ) such that (-A) 7/2 t; e L q (R N ). 
Proof. Using property (13. ip and Lemma I3.1[ we get 

(M"-\0(-A) 7/2 t; = / (-A)T/ 4 (Kf-\;)(-A)^ 4 t; 



th I ^MlEW^E^VE^)) 
Jr n+1 

^-K^l !/-|V(|E(. 

<r Jr* +1 



\q/2\\2 



> 4 (g - v f |(_Ar/>i*/ 2 | 2 



In the last step we get only inequality because the function | E(f)| 9//2 is not necessarily 
7-harmonic. □ 

With the same technique a generalization of (15. ip can be proved. 
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Lemma 5.2 Let < 7 < 2. T/jen 

(5.2) f ij(v)(-Ay /2 v > [ |(-A)^(tO| 2 
whenever ip' = (\I/') 2 . 

Proof. Use the extension method and the property (Vip(w), Vu>) = |V\I/(wj)| 2 . □ 

In order to prove our second inequality we need the well-known Hardy-Littlewood- 
Sobolev's inequality [25], [37], [SI]: for every w such that (-A)^/ 2 t; G L^E^), 1 < 
r < AT/7, < 7 < 2, it holds 

(5.3) |M| n <c(7V,r, 7 )||(-A)^;|| r , 

r i = n-~ r - Putting for instance r = 2, 7 = cr/2, we obtain the inclusion .H' tr / 2 (R Ar ) °->- 

1"-° (R^) whenever N > a. What happens for iV = 1 < a < 2? Or more generally, 
for r > iV/7? We answer this question in the next lemma. 

Lemma 5.3 (Nash-Gagliardo-Nirenberg type inequality) Let p > 1, r > 1, 

< 7 < min{iV, 2}. There is a constant C = C(p, r, 7, N) > such that for any 
v G L P (R N ) with (-A)~</ 2 v G L r (R N ) we have 

(5-4) \\v\\« 2 +1 <C\\(-Ay/ 2 v\\r\M a p , 

7 N(rp+r— p) p(r—l) 

where r 2 = — — \ , a = — -■ 

z r(N—-y) ' r 



Proof. We use (I5.ip . Estimate now the left hand side of this inequality using 
inequality (15. 3p . and the right hand side with Holder's inequality, to get (15 .4p . □ 

Notice that, for r = 2 and 7 = cr/2, this corresponds to using inequality (15. 3p for 
the space H a ^(R N ) instead of H°/ 2 (R N ), thus allowing all values of a G (0,2) even 
in the case N = 1. 

Inequalities of this kind are already available [8]. However, this particular one is, 
up to our knowledge, new. Let us explain in more detail the consequences of this 
inequality in relation to (I5.3p . 

Assume first that A^ > jr. Hardy-Littlewood-Sobolev's inequality (15. 3 p implies 
that if (— A) 7 / 2 t? G L r then v G L ri . Assuming also v G L p , then (15 .4p gives that 
v G L n fl L r ' 2 , which is always stronger that v G L Tl (1 LP . Both results coincide in 
the case p = r\ = r 2 . 

If on the contrary A^ < jr, we cannot apply Hardy-Littlewood-Sobolev, but (15. 4p 
gives that v G L p n L r \ 

We now consider the case of a bounded domain Q C R^. The characterization of 
the fractional Laplacian in terms of the extension to the half-cylinder allows us to 
repeat the proofs of Lemmas 15.11 and 15.21 in the case where the domain is bounded. 
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Lemma 5.4 Strook-Varopoulos' inequalities (15. ip and (15. 2p /ioW £r«e u>zt/i s?/6- 
stituted by Q C bounded. 

On the other hand, let w G i^Q^ 2 (fi). Consider its cx-extension w = E(i>), (cr = 27), 
and let m; be the extension of w by zero outside the half-cylinder. Then we have the 
estimate, see [10], H~9l. 



^ / y L ~ a \Vw\ 2 > S(a,N) / | Ti(w 



N+l 

+ 



2N 

JV-CT 



N-cr 
N 



i.e. 



(5.5) fi a / y^lVwI^^TV) / |v 

\</r 



2iV 

iV-CT 



N-cr 
N 



The left hand side equals || (— A)^/ 2 ^|||. That is, we have obtained inequality (15. 3 p in 
the case r = 2. 



From this point, we can repeat the proof of Lemma 15.31 which only uses the case 
just proved and Holder's inequality, thus obtaining inequality (15. 4p also for a bounded 
domain. 

More important is the following application. 

Lemma 5.5 (Sobolev type inequality) Let Q C M. N be a bounded domain, and 
let v be such that (-A)"</ 2 v E L r (tt), N > 1, < 7 < 2. Then we have 

(5.6) \\v\\ q <C(q,r,N n )\\(-Ay/ 2 v\\ r 

for every 1 < q < j^J if N > -yr, or for every q>l if 7 < N < 'jr . 

Proof. If N > 7r we just apply (15. 4p with p = , and apply Holder's inequality 
for the exponents 1 < q < j^f ■ If 7 < N < jr, and given any q > we 
apply (15.41) . this time with s = jj^ 9 < ^ — r histead of r, and p = q. We end again 
with Holder's inequality for the exponents 1 < q < □ 



6 Uniqueness 

In this section we prove the uniqueness parts of Theorems 12.11 and 12.31 

Notations. We will use the simplified notation u m for data of any sign, instead of 
the actual "odd power" |u| m_1 u. In the same way, w 1 ' m will stand for \w\™~ 1 w. In 
addition, || • || p will denote the norm || • ||jy>(Rjv) or || • ||lp(^)- 
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6.1 m > m*. Uniqueness of weak solutions 

Theorem 6.1 Let f G L 1 (R JV ) and m > m*. Problem (II. ip /ias at most one weak 
solution. 



Proof. We adapt the classical uniqueness proof due to Oleinik [M]. This will require 
u G L m+1 (M. N x (0, T)), which will be true if m > m*. To prove this we apply Holder's 
inequality twice, first in space and then in time, using inequality (15. 3ft . Assume first 
N > a. We have 



T 

JR N 



U 



m+1 



dxdt < 



u\ dx ) ( / |w| N ~ tT dx 



1-0 



dt 



P 



< CT max \\u(-,tW 

tG[0,T] 



2Nm 

\u\ N - a dx 



dt 



1-7 



where /3 



N(m~l)+a(m+l) 
N(2m-l)+a 



and 7 



N(m-l)+a 
N{2m-l)+a 



. Observe that m > m* implies 0, 7 G 



(0, 1). Applying now inequality fl5.3j) . we get 



it 



m+l 



Jk n 



JO 



K(-,t)||| CT/2 c2t 



1-7 



< C. 



In the case N — 1 and 1 < cr < 2 the computation is similar. But we have to use the 
Nash-Gagliardo-Nirenberg type inequality (I5.4[) instead to get the same conclusion. 
What we get in this case is 

7 / pT . p . 2 _ a \ 1-7 



T 

Jr n 



U 



m+1 dxdt < 



\u\ dx) 1 dt 



2m+l 



< CT 7 max ||u(-,t)||; 

te[0,T] 



<S+l-7 



|m| 2 " ct dx 
u m (-,t)\\l„ /2 dt 



dt 



1-7 

< C, 



where now 5 



<r(m+l)-l 



and 7 



m— 1+tr 



2m-l+cr ' 2m-l+cr' 

We now proceed with the core of the proof. Let u and u be two weak solutions to 
Problem (11. ip . We take the following function as test in the weak formulation 



<p(x,t) 



{u m -u m )(x,s) ds, 



< i < T, 



with </? = for t > T. We have 

(« - 5) (x, t) (w m - u m ) (x, t) dxdt 



Jm N 



T 



(-Ay /4 (u m -u m )(x,t) I (—A) a/4 (u m - u m )(x, s) dsdxdt = 0. 

'0 ^t^ 

Integration of the second term gives 

(m - 5) (x, t) (w m - w m ) (x, t) dxdt 



^1" 



+i / (f (-A) a /\u m -u m ){x,s)ds] dx = 0. 

2 Jrn \J J 
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Since both integrands are nonnegative, they must be identically zero. Therefore, 
u — u. □ 

Remark. The same proof works without any restriction on the exponent m provided 
u G L m+1 (R N x (r,T)). 

6.2 m < m*. Uniqueness of strong solutions 

Weak solutions satisfy the equation in (II .11) in the sense of distributions. Hence, if 
the left hand side is a function, the right hand side is also a function and the equation 
holds almost everywhere. This fact allows to prove several important properties, 
among them uniqueness for m < m*, and hence motivates the following definition. 

Definition 6.1 We say that a weak solution u to Problem (II. ip is a strong solution 
if d t u G L°°((t, oo) : L 1 (M JV )), for every r > 0. 

In the case of strong solutions the uniqueness result also provides a comparison 
principle. The following uniqueness proof is valid for all values of m > 0. 

Theorem 6.2 Let m > 0. If U\, u 2 , are strong solutions to Problem (II. ip with 
initial data f±, f 2 G L 1 (M Ar ), then, for every < t\ < t 2 it holds 

(6.1) / (ui - u 2 ) + (x,t 2 ) dx < / («i - u 2 )+{x,t 1 ) dx. 

Jr n Jr n 

Proof. Let p G C^R) n L°°(R) be such that p(s) = for s < 0, p'(s) > for s > 
and < p < 1, and let j be such that f = y/p 7 , j(0) = 0. We will choose p as an 
approximation to the sign function. 

Let us first assume that ti > 0. We subtract the equations satisfied by U\ and u 2 , 
multiply by the function <p = p(n™ — u™\ and integrate by parts to get 

1 1 d{ui ~ u2) P (uT - o = - r / (-a)-/ 2 « - <)p« - o- 

ti JR N 01 Jt! JWL N 

We now apply the generalized Strook-Varopoulos inequality (15. 2p . to get 

r d( Ul - U2 ) p{< _ uT) ^ _ c r* j |( _ Ar /4 i(< _ <)|2 < 

We end by letting p tend to the sign function. The case t\ = is obtained passing to 
the limit. □ 

For the problem posed in a bounded domain, the above proofs of uniqueness work 
without any change. 
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7 Existence for bounded initial data 




Crandall and Pierre developed in [22] an abstract approach to study evolution equa- 

tions of the form d t u + Aip(u) = when A is an m-accretive operator in L 1 and p is a 
monotone increasing real function. This allows to obtain a so-called mild solution us- 
ing Crandall-Liggett's Theorem. Our problem falls within this framework. However, 
such an abstract construction does not give enough information to prove that the 
mild solution is in fact a weak solution, in other words, to identify the solutions in a 
differential sense. We will use an alternative approach to construct the mild solution 
whose main advantage is precisely that it provides enough estimates to show that it 
is a weak solution, and later that it is strong. 

In order to develop the theory for Problem (11. ip . we will approximate the initial data 
/ by a sequence f n £ L 1 (M JV ) D L°°{R N ), and use a contraction property in order to 
pass to the limit. Hence, our first task is to obtain existence for integrable, bounded 
initial data. This is the goal of the present section 

We will construct solutions by means of Crandall-Liggett's Theorem [21], which is 
based on an implicit in time discretization. Hence, we will have to deal with the 
elliptic problem 



Equalities on M. have to be understood in the sense of traces. To show existence of 
a weak solution for this problem we approximate the domain by half-cylinders, 

B R x (0, oo), with zero data at the lateral boundary, dB R x (0, oo). We recall that 
in the case a = 1 a similar construction is performed in [3S], though there half-balls 
were used instead of half-cylinders. 

Now we have two choices: either we first pass to the limit in the discretization, to 
obtain a solution of the parabolic problem in the ball Br, and then pass to the limit 
R — > oo; or we first pass to the limit in R to obtain a solution of the elliptic problem 
in the whole space and then pass to limit in the discretization. We will follow both 
approaches (each of them has its own advantages) and will later prove that both of 
them produce the same solution. 

Instead of just considering problems in balls we will analyze the case of any bounded 
open domain Q C Mr, since it has independent interest. 

7.1 Problem in a bounded domain 

In order to check that the hypotheses of Crandall-Liggett's theorem hold, we have 



to prove existence of a weak solution w of (17. ip (defined in the standard way) and 
contractivity of the map g i— > w l ^ m {-, 0) in the norm of L^iVL) for the elliptic problem 
for all g £ L°°(fi). 



(7.1) 




in R^+\ 
on R N . 
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Theorem 7.1 Let Q C M. N be a bounded domain. For every g e L°°(Q) there exists 
a unique weak solution w £ Xq{Cq) to Problem (17. ip . It satisfies ||u>(-, 0)||oo < 
Moreover, there is a T -contraction property in L 1 : if w and w are the solutions 
corresponding to data g and g, then 

(7.2) I [w 1/m (x,0)-w l/m (x,0)] dx< [ [g{x)-g{x)] + dx. 



Proof. The existence of a weak solution, i.e., a function w £ Xq(Cq) satisfying 

(7.3) ^ [ y 1 ^(Vw,V<f) + [ u 1/n V- f 9<P = 0, 

Jc n Jn Jn 

v = Tr(w), for every test function ip, is obtained in a standard way by minimizing 
the functional 

+ f\v\^- fvg. 

z Jen m + 1 Jn Jn 



This functional is coercive in Xfi (Cq). Indeed, the first term ||iu||jf<r, and is using 
Holder's inequality, we have 

vg 

Now, the trace embedding (15. 5 p implies 



v. 



< \\v\\ ™ \\g\\ 2N < e\\v\\ 2 _2N_ + -||.gf|| 2 2jy . 

JV-ct N+cr jv_ ct 5 M+a 



(7.4) J{w)>C x \\wf x «-C 2 . 

For A r = l<cr<2we use inequality (15. 6ft instead. In fact, putting q = r = 2, 
7 = er/2, we get > C||^|| 2 - We obtain again (I7.4p . 

We now establish contractivity of solutions to Problem (17. ip in L l (Q). Let w and w 
be the solutions corresponding to data g and g. If we consider in the weak formulation 
the test function tp = p(w — w), where p is any smooth monotone approximation of 
the sign function, < p(s) < 1, p'(s) > 0, we get 

/jj<j I y 1 ~ a p'(w-w)\V(w-w)\ 2 + [ (w 1/m -w 1/m )p(w-w)- [ (g-g)p(w-w) = 0. 
Jc n Jn Jn 

Passing to the limit, we obtain 

f ( w Vm _ ~l/m )+ < h g -g) sign ( u; _ ffi) < f(g-g) + . 
Jn Jn Jn 

In particular, under the assumption g > we have w(-, 0) > 0. Standard comparison 
gives now w > in Cq. In the same way we can establish a contractivity property for 
subsolutions and supersolutions to the problem with nontrivial, on Q x (0, oo), bound- 
ary condition. We thus may take the constant function g = \\g\\oo as a supersolution, 
to get llu'^lloo < Halloo- We also deduce the estimate 

(7.5) ^ [ y l ~ a \Vw\ 2 < [ gw< C. 

Jc n Jn 

19 



□ 

We now construct a solution to the parabolic problem in the bounded domain. 

Theorem 7.2 Let fz C M. N bounded. For every f £ L°°(fl) there exists a weak 
solution (u,w) to Problem (13. 5p with u(-,t) £ L°°(fz) for every t > and w £ 
L 2 ([0, T]; Xq (C^)). Moreover, the following contractivity property holds: if (u,w), 
(u, w) are the constructed weak solutions corresponding to initial data f, f , then 



(7.6) / [u(x, t) - u(x, t)]+ dx < / [f(x) - f(x)]+ dx. 

Jn Jn 

In particular a comparison principle for constructed solutions is obtained. 

Proof. Crandall-Liggett's result only provides us in principle with an abstract type of 
solution called mild solution. However, we know that f !7.3 jl and (17. 5p hold, from where 
it is standard to show that the mild solution is in fact weak, see for example [33] . 
We recall the main idea: For each time T > we divide the time interval [0, T] in 
n subintervals. Letting e = T/n, we construct the pair function (u £ ,w £ ) piecewise 
constant in each interval (tk-i,tk\, where tk = he, k — 1, • ■ ■ , n, as the solutions to 
the discretized problems 

L a w £yk = in Cq, 

dw 



e.k 1/m 

t: = wJh — u E k-i on iZ, 



with u £ ^k-i = W£fc™i(-,0), u £) o = f. The mild solution is obtained by letting e — > 0. 
We still have to check that it is a weak solution. 

Crandall-Ligget's Theorem gives that u £ converges in /^(fz) to some function u £ 
C([0, oo) : L 1 (fi)). Moreover, ||w e ||oo < 11/llSo- Hence, w £ converges in the weak-* 
topology to some function w £ L°°(Cn x [0, oo]). On the other hand, multiplying the 
equation by w £j k, integrating by parts, and applying Young's inequality, we obtain 

- / \u £ (x,t)\ rn+1 dx + y, a [ [ y l - a \\7w £ (x,t)\ 2 dxdydt 
) Jn Jo Jc n 

<7-^TV f \m\ m+1 dx. 
(m + 1 



[m 



Passing to the limit, the following estimate is obtained for the weighted norm of |Viw|, 

fi a I [ y x - a \Vw(x,t)\ 2 dxdydt < 1 / \f(x)\ m+1 dx, 
Jo Jc n \ m + 1) Jn 

and therefore w £ L 2 ([0,T}; Xq(Cq)). Now, choosing appropriate test functions, as 
in [35], it follows that we can pass to the limit in the elliptic weak formulation to get 
the identity of the parabolic weak formulation. 

The contractivity obtained in each step is inherited in the limit. □ 
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7.2 The problem in the whole space 



Theorem 7.3 For every f E L 1 (IR Ar ) D L°°(R N ) there exists a weak solution (u,w) 
to Problem fl3"3|l . This solution satisfies u(-,t) E L 1 ^) n L°°(R N ) for every t > 0, 
and w E L 2 ([0,T] : X a (JBL+ +1 )) . Moreover, the following contractivity property holds: 
if (u,w), (u,w) are the constructed weak solutions corresponding to initial data f, f, 
then 

(7.7) / [u(x,t)-u(x,t)]+dx< [ [f(x)-f(x)] + dx. 

JR N JR N 

In particular a comparison principle for constructed solutions is obtained. 

Proof. Let us comment briefly the two constructions of the solution. 

Take to begin with as domain Q = Br, the ball of radius R centered at the origin, and 
let wn(fji) be the corresponding solution to problem (17. ip with datum fn = f- Xb r - 
Passing to the limit R — > oo we obtain a weak solution w to the elliptic problem in the 
upper half-space R+ +1 . Now we can follow the technique described above using the 
time discretization scheme thus obtaining a weak solution, whose trace on {y = 0} we 
call U = U(f). This is the mild solution produced by the Crandall-Liggett theorem 
and as such is unique. The contractivity property (17. 7p follows from (17. 6p . 

As to the second construction, we use the weak solution ur^r) to the parabolic 
problem posed in the ball Br, as obtained in Theorem 17.21 In the study of this limit 
we first treat the case where / > and fR approximates / from below. Then, the 
family of solutions ur^r) is monotone in R and also ur^r) < U(f) follows from 
simple comparison. In this way we ensure the existence of 

u(x,t) = lim UR(x,t) 

J?->oo 

It is easy to prove that u is also a weak solution with initial data / and u(f) < U(f). 
The equivalence of the two solutions depends on the already proved uniqueness result, 
see the remark after Theorem 16. 1[ 

In the general case where / changes sign, we use comparison with the construction 
for / + and /~ to show that the family Ur(/) is bounded uniformly and then use 
compactness to pass to the limit and obtain a weak solution. Again we end by 
checking that u — U. □ 

Remark. Since we have uniqueness, as a byproduct of the limit processes of con- 
struction we get the following estimates for weak solutions with initial data / E 
L 1 (R JV ) fl L°°(R N ): 

IK,t)||l<||/||l, |K,t)||oo< H/lloo. 

Remark. In the course of the proof we have obtained a unique weak solution to the 
nonlocal nonlinear elliptic problem 

(7.8) u+(-A) a/2 u m = g mQ 
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for every g £ L\n) n L°°(ft), both for = R N and for a bounded domain with 
homogeneous Dirichlet condition. This weak solution satisfies u £ L 1 (f2) D L°°(Vt). 



8 Existence for general data 

We prove here existence for data / £ L l (M. N ). The idea is to approximate the initial 
data by a sequence of bounded integrable functions and then pass to the limit in 
the approximate problems. The key tools needed to pass to the limit are the L 1 - 
contraction property and the smoothing effect. As a preliminary step we must show 
that the approximate solutions are strong. 

8.1 Strong solutions 

We prove here that the bounded weak solutions constructed in the previous section 
are actually strong solutions. We remark that the proof does not require boundedness 
of the solutions, but a control of the L 1 norm of the time-increment quotients. Hence 
the property will be true for the general solutions constructed next by approximation. 

We start by proving something weaker, namely that d t u is a Radon measure. 



Proposition 8.1 Let u be the weak solution constructed in Theorem \7.3[ Then dtu 
is a bounded Radon measure. 

Proof. If m = 1, a direct computation using the representation in terms of the kernel 
yields 



du , . 



All, 

l at 



If m ^ 1, following step by step the proof for the PME case, see [I] or [39], we get 
that the time-increment quotients of u are bounded in L 1 (M. N ), 



f8.2) 



u(x, t + h) — u(x, t) 



dx -W^Wt 11/1,1 + o(i: 



as h — > 0. Hence, the limit d t u must be a Radon measure. □ 

The next step is to show that the time derivative of a certain power of u is an Lf oc 
function. 

Lemma 8.1 The function z = satisfies d t z £ L^ oc (IR Ar x (0, oo)). 
Proof. If we could use d t w as test function, we would obtain 

-Hv f [ y 1 ~ e7 {V(dw/dt),Vw}dxdydt= [ [ (d t z) 2 dxdt, 
Jo Jr* +1 Jo Jm 
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from where we would get d t z G L 2 {JL N x [0, T]). Though d t w is not admissible as 
a test function, we can work with the Steklov averages to arrive to the same result, 
following [6]. For any g G L 1 1 oc (M Ar ) we define the average 

j rt+h 

9 h (x,t) = - J g(x,s)ds. 

We see that 

o h g{x,t + h) - g(x,t) 

dt9 ~ h 

almost everywhere. Since for our solution we have d t u h G L 1 {R N ), we can write the 
weak formulation of solution in the form 

d t u h (p dxds = —\x a I I y x ~ a (Vw h ,V(p) dxdyds. 

JR N JO Jr^ +1 

Let now choose the test function, tp = (d t w h , where ( = ((t) > 0. Then the above 
identity becomes 

I (d t u h d(u m ) h dxds =-n* f [ y 1 - a C(Vw h ,Vd t w h ) dxdyds 
o Jr n Jo Jr^ +1 

= f I y^C'lVw^ 2 dxdyds <C. 
2 Jo Jr^ 

We end by using the inequality (u m ) h u h > c[(u n! f 1 ) h ] 2 , see for instance [35], and 
passing to the limit h — > 0. □ 

We finally prove that u is an L\ oc function, and hence that dtu G L°°((r, oo) : L 1 (M iV )) 
for all r > 0. Therefore, a is a strong solution. 

Theorem 8.1 If u is a weak solution to Problem fll.ip such that dtu is a Radon 
measure, then u is a strong solution. 

Proof. The first step is to prove that dtu G L\ oc {M. N x (0, oo)). This follows imme- 
diately from Theorem 1.1 in [6], once we know that dt {u 1 ^ J G Li oc (M'' v x (0, oo)), 
see Lemma 18.11 Having proved that dtu is a function, estimate (I8.2p implies 

du 



while we have the estimate (18. ip for m = 1. □ 

We end this subsection with two more estimates that will be useful in the sequel. 

Multiplying the equation by u m and integrating in space and time, (recall that u is 
a strong solution), we obtain 
(8.3) t 

( [ \(-Ay/ 4 u m \ 2 dxds + ^— [ \u\ m+1 (x,t)dx = [ \f\ m+1 dx. 

Jo Jrn m + 1 J RN m + 1 J RN 
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Thus, we control the norm in L 2 ((0, oo) : H a / 2 (M. N )) of u m in terms of the initial 
data. 

Proposition 8.2 In the hypotheses of Theorems \2.1\ or Wlk iff 6 L m+1 (R. N ), then 
the solution to Problem (11. ip satisfies 

i 

m+l 



4) / Il« m (-,*)IIW*<— rrWmti- 

Jn m + l 



Another easy consequence of (18.31) is that the norm ||it(-,£)|| m+1 is nonincreasing in 
time. In fact, this is true for all LP norms. 

Proposition 8.3 In the hypotheses of Theorems 12.11 or I2.3|. any LP norm, 1 < p < 
oo, of the solution to Problem (11.11) is nonincreasing in time. 

Proof. We multiply the equation by \u\ p ~ 2 u with p > 1, and integrate in ~R N . Using 
Strook-Varopoulos inequality (15.11) . we get 

d 



dt 



u\ p {x,t)dx =-p {-&y /2 {\u\ m - 1 u)\u\ p - 2 u 



< -C 



, CT /4i , ™+P-l 2 



-A) CT />| 



< 0. 



The limit cases p — 1 and p = oo were obtained from the elliptic approximation. □ 

Remark. The previous result can be easily generalized substituting the power \u\ p by 
any nonnegative convex function ip(u), using (15. 2p . Thus we obtain 

^[ iP(u)(x,t)dx < - [ |(-A) CT/4 ^(n)| 2 < 0, 

dt J JJJV J R N 

where \1/(m) = Jjj"' ^/ms m ~ 1 ^"(s) ds. 



8.2 Smoothing effect 

As a first step we obtain a bound of the L°° norm in terms of the LP norm of the initial 
datum for every p > 1, with the additional condition p > p*(m) = (1 — m)N/a if < 
m < m*. The important observation, that will be used in the next subsection when 
passing to the limit for general data, is that the estimates do no depend qualitatively 
on the L°° norm of the initial value. 



Theorem 8.2 Let < a < 2, m > 0, and take p > max{l, (1 — m)N/a}. Then for 
every f e L 1 (R N ) n L°°(R N ), the solution to Problem §L1§ satisfies 

(8.5) sup \u(x,t)\<Ct-*»\\f\$ 

x£R N 

with 7 P = (m — 1 + crp/N)^ 1 and 5 P = apj p /N, the constant C depending on m, p, a 
and N . 
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Proof. We use a classical parabolic Moser iterative technique. 

Let t > be fixed, and consider the sequence of times t k = (1 — 2~ fc )l We multiply 
the equation in (11. lj) by |-u| Pfc ~ 2, u, > Po > 1) and integrate in WL N x [t fc ,t fc+1 ]. As in 
the proof of Proposition 18. 3[ using (15. ip and the decay of the LP norms we get 



k 

k+1 



Pk II "A ' / llpfc 

The constant c4 depends on p ( as we h as on m an d AT, but not on a). We now use 
the Nash-Gagliardo-Nirenberg type inequality (15. 4 p with r = p k + m — 1 an use again 
the decay of the LP norms to obtain 

' fc+ ^l<,r)||^||(-Ar/ 4 |n| Qd F i (-,r)||^r > C H- : r) f^XU <*r 

t k Jt k 2N-<7 



> C2- k t\\u(;t k+1 ^ 2pk+m - ] 



+ m— 1) * 
2JV-ct 



Summarizing, we have 



\\u(-,t k+1 )\\ Pk+1 < (2 k d , k t- i r^\\u(- J t k ) 



s pk 



Wk J 



where = s(p fe + ^^), s = > 1. 

First of all we observe that taking as starting exponent po = p > (1 — m)N/cr (and 
p > 1) it is easy to obtain the value of the sequence of exponents, 

p k = A(s k - 1) +p, A = p-± '— > 0. 

a 

In particular we get p k +i > p k , with Um^ooPk = oo. Observe also that min{l,m} < 
— — - < max{l,m|. This implies that the coefficient in the above estimate can be 

Pk+rn— 1 — l J J 1 



bounded by c Pfc + x , for some c = c(m,p, N, a). Now, if we denote U k = \\u(-,t k )\\ Pk , we 
have 

k s J*h- 

Pk+i 



U k+1 < c^+i t 2p k+i u k 
This implies U k < c ak t~ l3k UQ k , with the exponents 

1 STn i i N(N-a) 1 A JV s*p p 

ai = _p- j)s ^__, A = _g^_ 4 = -^ 

We conclude that 



|u(-,*)lloo= lim u k <cr%u* A =cr(™-"»***\\f\\$ 



pa 

(m — TJW+po 
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□ 



Remark. When a < N, which is always the case if N > 2, we may use the Hardy- 
Littlewood-Sobolev's inequality ( 15. 3 p instead of the Nash-Gagliardo-Nirenberg type 
inequality ( 15. Ah to arrive at the same result. 

The constant in the previous calculations blows up both as p — > p*(m), < m < m*, 
or as p — > 1 + in the case m > m*. Nevertheless, in this last case an iterative 
interpolation argument allows to obtain the desired L l -L°° smoothing effect. 



Corollary 8.1 Let < a < 2, m > m*. Then for every f e L 1 ^) n L°°(R N ) , the 
solution to Problem ( 13. 5 P satisfies 

(8.6) sup \u(x,t)\ <Ct^\\f\\{ 

with 7 = 71 = (m — 1 + cr/iV)" 1 and 5 = 5\ = a^/N, the constant C depending on 
m, N and a. 



Proof. Putting = 2 k t, estimate ( 18 .5p with p = 2 for instance (for which it is 
valid if m > m*), applied in the interval [t\, to] gives 

2(T72 CT 72 °"72 

\\u(-Moo<c(t/2)-'*\\u(-,T 1 )\\^ <c(t/2)-^||n(-,r 1 )|| 1 ~||n(-,r 1 )H^ 



oc 



We now apply the same estimate in the interval [t2, T\], thus getting 



°"72 

7 72 / 20-72 \ 



||u(.,*)lloo < C {t/2)-^\H;ri)\\7~ C(t/4)-^\\u(;r 2 y^ 



12 

Iterating this calculation in [Tfc,Tfc_i], using Proposition 18.31 we obtain 

n«(-,t)iioo < c a *2 b H- d *\\ U (-Mn<-,n)\\i k - 

Using the fact that m > m* implies ^ = ^ m _-^ N+2(T < 1, we see that the exponents 
satisfy, in the limit k — > oo, 

- H 2<7 V (y»~ l)iV + 2a _ 7iff 
ak ~^N) ( m -l)N + a N + ' 

m-l)7V + 2<7 



fc = 5>«-H>(£) : 

j'=o 

4 = 720fc -> 7 ryr— = 71 

efc = a fc - 1 -)> 



((m- l)iV + a) 



2 



A = 2(^ 



iV 7 



k 

->■ 0. 
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□ 

Remark. Since the LP norm is nonincreasing, an U decay is obtained again by 
interpolation, for any r > p > 1, 



(r — p)N 



\U\ 



<^ q£ r((m—l)N+pa) 



p((m—l)N+r(7) 
r((m — l)N+pa) 

V 



8.3 Passing to the limit. Existence of strong solutions 

We first prove the existence of a strong solution for the case of initial data with 
improved regularity 

Theorem 8.3 Let a G (0, 2) and m > 0. Then for every f G L 1 (M 7V ) if m > m* or 
f G L 1 (M 7V ) (1L P (M. N ) with p > p*(m) = (1 —m)N/cr if m < m„ there exists a strong 
solution to Problem (II. ip . 

Proof. Let {fk} C L 1 (M iV ) D L 00 (R JV ) be a sequence of functions converging to / 
in L 1 , and let {uk} be the sequence of the corresponding solutions. Thanks to the 
/^-contraction property, we know that Uk{-,t) — > u(-,t) in L 1 (IR Ar ) for all t > for 
some function u. Moreover, nonlinear Semigroup Theory guarantees that Uk — > u in 
C([0,oo) : L 1 (R iV )) [20], [23]. 

Consider a fixed time r > 0. Using the smoothing effect and the estimate (18. 4ft . we 
have u k G L°°(R N x [r,oo)) and u 1 ^ G L 2 ((r,oo) : H a/2 (R N )), both uniformly in fc. 
Thus the limit u is a weak solution to Problem (II .ip for every t > r. We now want 
to go down to t = 0. This follows from the /^-contraction and the /^-continuity. In 
fact 



\u(x,t) — f(x)\ dx < / \u(x, t) — Uk(x, t)\ dx + \ \uk(x, t) — fk(x)\ dx 
+ [ \f k (x)-f(x)\dx. 



The fact that u is a strong solution is now an immediate consequence of Theo- 
rem O □ 

As a byproduct, and using the uniqueness results, Theorems 16.11 and 16 . 2 1 we obtain 
the /^-ordered contraction for the limit solution: given u and u two weak solutions 
to Problem (II. ip . then for every < t\ < t 2 we have 

(8.7) / [u(x, t 2 ) — u(x, t2)]+ dx < / [u(x, ti) — u(x, dx. 

Jr n Jr n 

We now consider the case of general data, / G L 1 (M Ar ). We only need to look at the 
case < m < m*, in view of the previous theorem. We prove here that the unique 
mild solution is in fact a very weak solution. 
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Theorem 8.4 Let a G (0,2) and m < m*. Then for every f G L 1 (R Ar ) there is a 
unique mild solution to Problem (II. ip which is moreover a very weak solution. 

Proof. As in the proof of the previous theorem, we approximate the initial data 
by a sequence {fk} of bounded, integrable initial data. The corresponding strong 
solutions Uk converge in C([0, oo) : L 1 {R N )) to a certain function u which, being the 
limit of mild solutions, is also a mild solution according to the theory, [21]. Moreover, 
integrating by parts in space, we see that the approximate solutions Uk are very weak 
solutions, namely 

(8.8) rf u^-Tf <vAr/v = o 

Jo Jw N at Jo Jr n 

for every smooth and compactly supported test function (p. On the other hand, 
||ujTlli/m = IK!!™ < \\f\\T- Hence uf -± u m in L 1 ^ m (R N ). Since (-A) CT /V belongs 
to the dual space LV^-^iR"), we conclude that u is a very weak solution. □ 

The passage to the limit in the case where the spatial domain is bounded is similar. 



9 Further qualitative properties of the solutions 

We prove in this section some important properties that our solutions have. Through- 
out this section u is the strong solution to Problem (II. ip corresponding to an initial 
value / satisfying the hypotheses of Theorem 18.31 



9.1 Positivity and regularity 

We start with this observation: if u > is a classical solution and u(xo,t) = for 
some Xq and t, then formula (11.31) gives (— A) a ^ 2 u m (x , t) < 0, unless u(-,t) = 0, and 
hence dtu(xo, t) > 0. Therefore, we expect solutions with nonnegative data to become 
positive immediately, and to stay positive unless they vanish. However, solutions are 
not known to be classical. Hence, we will use a different argument, which involves 
the extension Problem (13.51) . 

The first ingredient in our proof is a control of the decay of the solutions, that has 
an independent interest. This control is based in an argument of Alexandrov's type, 
cf. [2], which is a bit delicate in this case, since the function u is not yet known 
to be continuous. 

Proposition 9.1 Assume f has compact support. For any bounded measurable set 
M with \M\ > there is a large enough radius R*, depending only on the support of 
the initial data and on M , such that 

u(x,t) < sup u(z,t) for a.e. \x\ > i?*, t > 0. 

zeM 
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Proof. Let supp(/) C B(0,R), M C B(0,R'). Thanks to Lebesgue's density 
Theorem, we know that there is a point X\ £ M such that for all 5 > there exists a 
radius such that for all cubes Q C B(xi,rg) we have 

I"n0| 



IQI 

Let us now consider the cube Q* = Q(x*, 2), with x* far away from the origin to be 
chosen later. We can cover Q* (except a subset of zero measure) with a finite number 
of disjoint cubes, small enough such that they are reflections of cubes centered at x\ 
and contained in B(xi,r$). Hence, an argument of Alexandrov's type (which can be 
done if \x\* is large enough) shows that 



\{u(-,t) < swp x&M u(x,t)}nQ* 

\Q*\ 



>l-5. 



More precisely, let Q be any of the cubes covering Q*, and let x be its center. We 
reflect around the hyperplane in Q, ir = (x — Xi) ■ (x — {x\ + x)/2) = 0, getting 
a cube Q' C B(x\,rs). It is clear that if \x\ is large (depending on R and R'), 
then the hyperplane ir divides the half-space Q in two parts, Q = Hi U H2 with 
S(0, R) x [0, 00) c H±, Q x [0, 00) C Hi- In this way, by the comparison principle, we 
obtain that the function z(x, y, t) = w(x, y, t)—w(x+xi — x, y, t) satisfies z(x, y, t) > 
almost everywhere in Hi, t > 0. 

Therefore, 

\{u(-,t) < supu(x,t)}nQ| > \{u(-,t) < sup u(x,t)}nQ'\ > (l-5)\Q'\ = (1-5)\Q\. 
Summing up, 

\{u{-,t) < supu(x,t)}nQ*\ > (1-5)\Q*\. 

We end by letting 5 — > to get that u(-,t) < swp xeM u(x,t) a.e. in Q* . □ 

The second ingredient in the proof of positivity is the following estimate on the time 
derivative for nonnegative solutions, arising from the homogeneity of the parabolic 
operator, see for instance [5], 

. du 

(9.1) (m - l)t— + u > 0. 

If m = 1 this formula is empty. However, in that case we have the representation 
formula (jl.5p . in terms of the fundamental solution K a , from which it is easy to derive 
the estimate 

atd t u + Nu > 0. 

It is sharp: equality holds for the fundamental solution at the origin. 
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Theorem 9.1 Let f > 0. Then for every t > 0: either «(•,£) = or mix u(x,t) > 
°k > for every compact K C M. N . 

Proof. By comparison, we only need to consider compactly supported initial data. 
Let then the support of / be contained in the ball B(0,R). Let K C M. N be any 
compact set, and assume that there exists some t > such that mix u(x,t ) = 0. 
Then there is a measurable set M e C iT, |M E | > 0, such that u(-,to) < e in M e . By 
Proposition 19.11 there exists some ball B(x*, 1) where u(-,to) < e almost everywhere. 
Letting e-^flwe get u(-,t ) = almost everywhere in B(x*, 1). Then we use (19.11) 
to prove that u is zero in that ball for an interval of times /. In fact, if m > 1 we 
obtain from (19.11) 

u(-M)>u(-M){fj 

a.e. in M N , ti > t\ > 0. Therefore we deduce u(-,t) = a.e. in B(x*,l) for every 
< t < t , that is, we may take I = [0, to]. In the case m < 1 formula (19. ip gives the 
same property in I = [t , oo). The case m = 1 follows from (II. 5p . 

We now observe that the integral definition of solution implies that for every test 
function ip that vanishes on dB(x*, 1) x (0, oo), we have 

0= II u^- dxdt = \i a I I I y l ' a {Ww,W^) dxdydt 

JiJb{x*,1) <Jt J I J0 JB{x*,l) 

dw , , 
— — ip dxdt. 

i jb(x*,i) vy a 

This gives d y vw(-, 0, t) — a.e. in x G B(x*, 1) for all t £ /. Now, for each fixed t £ / 
we extend w in a even way in the y variable to obtain a solution of the elliptic equation 
L a w = in a ball i?((x*,0), 1) C IR^ 4 " 1 . Since A(sc, y) = \y\ l ~ a is an A 2 -weight, we 
can apply a half Harnack's inequality (Theorem 2.3.1 in [23]), to obtain 

R /?L^ m ( M) - i n L M w &^ - c M-,t)\\ L 2 mx , m)) . 

B(x*,l/2) B((x* ,0),l/2) 

If u(-,t) ^ in M^, then w(-,t) > in f2, and this results in a contradiction. □ 

As a consequence of positivity, we can extend to the whole range m > the conti- 
nuity of solutions that was obtained in [3] for m > 1. 

Theorem 9.2 issume u > and u(-,T) 0. T/ien u £ ^(M" x (0,T)) for some 
< a < 1. 



Proof. The above-mentioned regularity result of [3] applies for bounded solutions to 
the equation 

+ (-Ay/ 2 v = 



dt 

in some ball B C and t > 0, with a nondegeneracy condition on the constitutive 
monotone function (5. This condition is fulfilled once we know that in any given ball 
the solution is essentially bounded below away from zero. On the other hand the 
solution is bounded for every positive time, thanks to the smoothing effect. □ 
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9.2 Conservation of mass 
Theorem 9.3 Let m > m*. Then for every t > we have 

u(x, t) dx = I f{x) dx. 



Proof. Thanks to the /^-contraction property, it is enough to consider the case of 
bounded initial data. The proof will follow different arguments in the cases m > m* 
and m = m*, and even in this latter case we have to distinguish between dimensions 
N > 2 and N = 1. 

Case m > m*. We can adapt the technique that was used to prove the property in 
the case a = 2 to deal with the nonlocal operator. It works as follows: we take a 
nonnegative non- increasing cut-off function ip(s) such that ip(s) = 1 for < s < 1, 
tf)(s) = for s > 2, and define ^>r{x) = if)(\x\/R). Multiplying the equation by ip R 
and integrating by parts, we obtain, for every t > 0, 

(9.2) ^ / WPR = ~ [ u m (-A) CT /V«. 



dt „ 

The radial cut-off function tpn has the scaling property 
(9.3) (-Ay/ 2 M%) = R- a (-AY /2 Mx/R). 

In addition, (— A)""/ 2 ^ E L 1 (K N ) fl L°°(M. N ). Both properties are straightforward 
using representation (jl.3p . 

Then, if we apply Holder's inequality with p = max{l, 1/m} to the right-hand side 
of (19. 2p . and use the above property, we get 



d 




dtj 





< !l/lir 1/p ll/lli /p ll(-A) CT/2 ^llp/(,-i) 

We conclude since the exponent of R is negative precisely for m > m*, and thus 



d f d f 

— I u = lim — / u<pn = 0. 

dt ./TD.JV R— dt J TON 



Case m = m*, N > 2. This case is much more difficult. The idea is to study 
separately the behaviour of the mass in a bounded set and close to the infinity, and 
also to decompose the fractional Laplacian into two operators. First, for every given 
5 > 0, we put u — «i + u 2 , where 

Ul = u ■ X{|x|<ij }, / < <5- 

Jr n 

Observe that u m = u™ + u™. Now express these functions in the following form 

u ™ = (-A) 7/2 z, u™ = ew + (-Ay /2 w, 
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where 7 = 2 — a > 0, e > 0. Then our equation becomes 

^ = Az + Aw - e(-Ay /2 w. 

We have introduced the e-regularization in the definition of w since, when applying 
inequality ( 15. 3p . it gives no information in the critical case m = m* if e — 0. 

As before, multiplying by the test function tpp>, and integrating by parts, we have 

-I / uip R <[ \zAip R \+ f \wA m \+e [ \(-A) l/2 w (-A^V-rI 
at Jr n Jm N Jr n Jm N 

= h + I 2 + h ■ 

We estimate each integral using Holder's inequality, the properties of (pn and some 
estimates on z and w. 

Estimate of h. We have G L r (R N ) for every 1 < r < oo, with 

, , (l-mr) iV(l-mr) 

IKIIr<IKIiri{N<^}l - <ci2o r ■ 

Using inequality (EJ) we have that z = (-A)-T/ 2 (n^) G L^fR^), g = — , for 

N — 7r 

every 1 < r < N/j, with ||z|| f/ < c||u™|| r . Then for such values of r, and using 
Holder's inequality and (I9.3p . we have 

JV(l-mr) N(l— mr\ 

\h\ < \\z\\ q \\&ip R \\m-i) ^ cR o r R- 2+N{q - 1)/q = c(R /R)^- . 
The exponent is positive if we take 1 < r < 1/m. 

Estimate of I<i. We have here u™ G L r {R N ) for every 1/m < r < oo, with 

II „ .m II II , II m ^ „cm 

||«2 II 1/m = HHIl ^ C ^ • 



If we multiply now the equation satisfied by w by u; 1 /" 1 1 ) and integrate in WL N , we 
get 

e [ w 1/m + [ w l/m ~ 1 {-Af /2 w = ! u™w 1/m -\ 
Jr n Jr n Jr n 

The second term is nonnegative by (15.11) . Thus, by Holder's inequality we get 

e|M|i/ m < IK1li /m |k 1/m |ll7 ( ?„ m) , 

i.e., e||iw||i/ m < H^Hi/m < c5 m . This implies 

|/ 2 | < |h||i/m||A^||i /( i- ro) < c^e- 1 ^- 2 . 

Estimate of J 3 . Since from the previous calculations we have ||(— A) 7 / 2 it;||i/ m < 
2||^||i /m < c5 m , we get 

\h\ < £||(-A)^|| 1/m ||(-A) CT -VB||i/(i-m) < c5 m eR 2 ^. 
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Summing up, we have obtained 



d I 




dt Jj 


IN 



jV(l-mr) 



< r(/,' () //.')- — + c5 m e~ l R a ~ 2 + c5 m eR 2 - a . 



We now choose e = R a 2 , and make first R — > oo and then 5 — > to conclude. 

Case m = m*, N = 1. Observe that m* > implies a < 1. Here we consider the 
same functions z and w as before, but with 7 = 1 — er > 0. The equation becomes in 
this case 

<*L = -(-A) 1/2 z - (-A) 1 /^ _ £ (-A) ff /V 

From here on, the calculations are exactly the same as the ones for the case m = m*, 
N > 2. □ 



9.3 Extinction 

The condition m > m* to have mass conservation is not technical, as shown by the 
next result on extinction in finite time, which extends the result by Benilan and 
Crandall for the standard differential case a = 2, see [5]. 

Theorem 9.4 Let < a < min{2, N} and < m < (N - a)/N. Then, if f 6 

L (l-m)N/a^N^ there ^ fl fi mte Ume T > q s ^ ^ U ( X ^T) = O.C 



Proof. From the proof of Proposition 

d 
It 

Using inequality H5 .31) we obtain 



< -C 



we have 
(_A)-/ 4 



m+p— 1 
U\ 5 



\u\ p dx + C 



AT-CT 

(y+m-l)AT \ -jf- 

ml cb I < 0. 



If we now choose p = (1 — m)N/a (this is where the restriction on m comes, since 
p has to be bigger than one), we get that the function J(t) = \\u(-,t)\\ p satisfies the 
differential inequality 

J 1 it) + C [t) < 0. 
This implies extinction in finite time provided J(0) is finite. □ 

Remark. In the special case m = (N — cr)/(N + a), N > a, there exists an explicit 
family of solutions in separated variables with extinction in finite time, 



u(x, t) = b Nj(7:C (T — t) 2jv [ c + \x — a\ \ 2 



T, c > 0, a G 



The spatial part satisfies the elliptic fractional equation (— A)°~/ 2 Lp m = see for 
instance [17] . 
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9.4 The problem in a bounded domain 

It is easy to see, following the ideas of the proofs for the case where the domain is 
the whole space, that the smoothing effect is true for the solutions to the problem 
posed in a bounded domain Q C M N (though the obtained decay rate will not be 
optimal for the problem in a bounded domain). The decay of the solution to yields 
the decay of the mass. Even more, solutions become extinct in a finite time for very 
< m < 1. 



Proposition 9.2 Let < a < 2 and < m < 1. Then, if f G L p (f2) for some 
p > rnax{l, (1 — m)N/o~}, there is a finite time T > such that u(x,T) = a.e. in 

n. 



Proof. With the same calculations used for the whole space, using here inequality 
(I5.4p . we obtain 

j t J \u\ p dx + c(J \u\ p dxY<0, 

where 7 = (p + m — l)/p G (0, 1). □ 

Retention. For m > 1 it is easy to see that nonnegative solutions do not become 
extinct in finite time, not only in the case of the Cauchy Problem posed in M. N , but 
also for the problem posed in a bounded domain with zero Dirichlet data. Indeed, 
if m > 1 nonnegative solutions satisfy estimate ( 19. ip . Hence, the function t m ~ 1 u is 
nondecreasing. This implies a retention property: if a solution is positive at some 
point at some time, it will remain positive at that point for any later times. 

The above retention property for nonnegative solutions is also true for the case 
m — 1, though a different proof is needed. Let O C Q be an open set in which 
u(-,ti) > (recall that u is continuous for t > 0). Let y? = (fo,i be the normalized 
eigenfunction corresponding to the first eigenvalue fi = fio,i of (— A) CT//2 in O with 
homogeneous Dirichlet boundary condition. We take e > small enough so that 
u(-,ti) > Eip a.e. in O. The solution to the fractional heat equation in O, t > t\, with 
initial data eip and homogeneous Dirichlet boundary condition is given by 

v(x,t) =ee-^ t - tl) tp(x). 

We claim that u > v for all x G O, t > t\, from where the retention property follows, 
since (f > in O. 

The claim is proved using a comparison argument in O. Comparison is in principle 
not at all obvious, since the fractional Laplacian operator changes with the domain. 
However, we have 

(-A) a J 2 u > (-A) a J 2 u inOcO. 

This is proved easily for any nonnegative u in the right spaces using the extensions of 
u to the corresponding half- cylinders. Hence d t u + (— A)^ 2 u m > 0, from where the 
claim follows. 
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An analogous argument can be performed for nontrivial nonpositive solutions to 
prove that they never become extinct if m > 1 . 

Remark. In the case m = 1 non-extinction is true for any nontrivial initial data 
in Indeed, the solution becomes bounded immediately after t = 0, and re- 

mains nontrivial for a while (recall that it is continuous in In particular, 

n(-,r) G L 2 (Q) and is nontrivial for some small r. Since u is smooth for posi- 
tive times (in fact C°°), the solution can be expanded as a series of eigenf unctions, 
«(',*) = EZi^e-^- T) ip k . Hence \\u(;t)\\ L 2 {a) = ^Zi \a k \ 2 e~ 2 ^-^ for all t > r. 
Notice that this quantity is positive for all t>r, since u(-,t) was nontrivial. 

10 Continuous dependence 

In this section we prove the continuous dependence part of Theorem 12.21 The result 
we obtain below includes the limit case m = m* for N > 2. The case o — > 2 is a bit 
different from the rest, and is dealt with separately. 

10.1 0<a<2 

We consider for N > 2 the region 

D = {(m, a) : < a < 2, m > (N - a)/N}, 
while for N = 1, 2 we take 

D = {(m, a) : < a < 2, m > {N - a)+/N}. 

Theorem 10.1 I7ie map S : D x L l (R N ) C([0,oo) : L 1 (M Ar )) zs continuous in 
all the arguments (m,a, /). 



This will follow from a result of nonlinear Semigroup Theory which states that if each 
of A n , n = 1, 2, . . . , oo is an m-accretive operator in a Banach space X, f n G D(A n ) 
and u n is the solution of 

—j^ + A n u n = 0, w n (0) = f n , 

then A n — > and f n — > foo imply u n — > Mqo in C([0, oo) : Af), where A n — )• A is 
understood as 

lim (/ + A^g = (I + Aoo)" 1 ^ for all g 6 A". 

See, e.g, [20], [23] for statements and references. Hence, the theorem will be a corollary 
of the convergence of {I + A mn ^ n )-\ where A m ^(u) = (-A) a / 2 u m . Thanks to the 
contractivity in L 1 (M 7V ) of the elliptic problems under consideration, it is enough to 
prove this convergence for functions g which are also bounded. 
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Proposition 10.1 Let {(m n , a n )}^L v (m n ,a n ) G D, be such that m n — » m and 
a n -> a as n ->■ oo, (m, a 7 ) G D. Then, for all g G L 1 (R N ) n L 00 ^) 

lim (J + A m ^J^g = (/ + A m ^)~ l g. 

n—toc 

Proof. Step 1. Let u n = (I + A^^) -1 ^ i.e., w n is the unique solution to the 
equation 

(10.1) u n + (-Ay^ 2 u^=g, 

see (]7.8p . The /^-contraction estimate (17. 2p . which is also valid for the case where 
T = R JV , implies the bounds 

\\u n \\i < \\g\\i, 

\\u n ~ ThUnWx < \\g - T h g\\\ , 

for each h G M. N , where (Thv)(x) = v(x + h). This is enough, thanks to Frechet- 
Kolmogorov's compactness criterium, to prove that {u n } is precompact in L l (K) for 
each compact set K C M, N . 

Step 2. To extend compactness to the whole M. N we need to control uniformly 
the tails of the solutions at infinity. More precisely, we need to prove that, given 
e > 0, there exists some R > such that ||wn||z,i(RJv\B H (o)) < £ - This will follow from 
a computation which is very similar to the one in the proof of Theorem 19.31 but now 
taking as test function 1 — ip^ instead of ip^. The problem is that this convergence 
fails to be uniform when m approaches m*. In order to reach m* we use a different 
argument, which unfortunately only works when N > 2. 

We first perform a reduction in order to consider only a characteristic function as 
initial value, a technique borrowed from [5]. Given e > 0, there exist M, ro > such 
that h = Mx{\ x \<ro} satisfies \\(g — h)+\\i < e. Notice that g < g + < h + (g — h) + . 
Hence, using the /^-contraction property (I7.2p . we have, denoting by Uh the solution 
corresponding to the initial data h, 



(Ug) < U g + < / U h + / \U h+{ g_ h)+ ~Uf\ 

x\>R J\x\>R J\x\>R J\x\>R 



< / u h + \\(g - h)+ ||i < / u h + e. 

J\x\>R J\x\>R 

Hence, it is enough to consider the special case where g is a regular, radially symmet- 
ric, approximation of Mx{\ x \<r }, with support contained in {|x| < r + 1}. Notice 
that the same property of radial symmetry is true for the solution for any later time. 

Step 3. We now write equation (110. ip as 



u n - Az n = g, z n = (-A)-( 2 -^/ 2 < 



v— (2— ct„)/2 ,m„ 

whose weak formulation is 

(10.2) / u n (p+ I (Vz n ,Vif)= I gtp 

Jr n 
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for every test function (p. For radial solutions this equation reads (abusing notation) 
r^V - (r^z'J = r N - x g (=0 for r > r + 1). 

Since u™ n G L P (M. ) for all max{l,l/m n } = p n < p < oo, and since N > 2 
implies m n > (2 — a n )/N, and thus p n < N/(N — a n ), applying ( 15. 3D we get that 
2 n = (-A)^ 2-47 ")/ 2 ^" G //(R*) for every N/(N - 2) < r < oo. On the other hand 
z n G C ,1,7 (M 7V ) for some 7 G (0, 1) depending on a = max n {a n } < 2 [36]. We even 
have, as in [11], the estimate H^nllc 1 ^ < c independent on n. On the other hand, the 
quantity $ n = r N ~ x z' n will be seen as a flux, and indeed 



00 

N-l 



r N u n (r) dr = \$ n (R) - ® n (oo)\ for every R > r + 1, 



which implies that the uniform control of the tails is equivalent to the uniform control 
of $ n (r) for large r. It is clear that there exists a limit flux, $ n (oo) = lim^oo <& n (r), 
since it is monotone increasing, 3>' n (r) = r N ~ l u n > for r > ro + 1. Notice also that 
integrating the equation, we get 



OO /"OO 

N-l.. f„\ A„. / „N-1. 



$ n (oo) - $ n (0) = / r N ^u n (r)dr- r N ~ l g(r) dr < 0. 



Since $ n (0) = 0, we get $ n (oo) < 0. Moreover, $„,(r) < for r > r + 1. 

Let us prove that $ n (oo) = 0. If there exist constants C, r\ > such that $ n (r) < 
— C for r > ri, after an integration we get 

^n(^) > Cr 2 ~ N for r > ro + 1. 

This is a contradiction with the property z n G L Ar// ( 7V_2 )(R 7V ). Thus $ n (oo) = 0. 

Step 4. We now have to estimate this limit more carefully. We claim that $ n (r) 
is small for large r uniformly in n, which will give the desired uniform control of the 
tails. 

Hardy-Littlewood-Sobolev's inequality ( 15. 3 p gives a control of the norm ||^ n || n in 
terms of the norm ||w™ n || Pn . But this norm can be estimated easily in terms of the 
norms \\g\\i and \\g\\oo an d the value m = max n {m n }, with constants independent on 
n. Therefore 

poo jy 

(10.3) / r N ~ l z%- 2 (r) dr < C 

Jo 

for every n > 1. 

N 

We now adapt some ideas from [5] . For e > given, if we take R(s) = r e 2Ce N ~ 2 , 
we have the estimate 

*m N 

r N ~\er 2 ~ N )—-dr>2C. 

'r 
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Integrating the inequality 

JV-l 



zM<(^) <(*), r + l<r<s, 



in s for s G [r, t], and taking t > 2r, we get, since z n > 0, 

\<$> n {t)\ = \t N ~ 1 z' n {t)\<cr N - 2 z n {r) ] 
i.e., \Q n {2R(e))\r 2 ~ N < cz n (r) whenever R(e) > r. Putting all together we get 



R(e) „ r R(e 



N 



cl r N - l (\§ n (2R(e))\r 2 - N )—idr < r N ~ l z^ 2 {r) dr 

r R(e) 

<C< r N -\er 2 - N )^dr, 

Jro 

which implies 

\$n{R)\ < \$„{2R(e))\ < ce, for every R > 2R(e), n>l. 
This ends the uniform control of the tails. 

Step 5. Summing up, we have obtained that along some subsequence, which we 
also call {(m n ,a n )}, the following convergence holds 



— > u* in L 



for some function it*. What is left is the identification of the limit, that is u* = m^j. 

The convergence u n — > u* in L 1 (WL N ) implies the convergence u™ n — > u™ in 
L r (R N ), r = max{l, l/min„{m n } }. Therefore, using inequality (15.30 . z n — > z = 
(_A)-( 2 ~-)/ 2 uf in Lt£*(R n ). 

We now take z n as test function in (110.21) . use that u n > 0, and apply Holder's 
inequality with exponents p = N/(N — 2),q = N/2, and obtain, thanks to (110.30 . a 
uniform control of the gradients of z n , 

|VzJ 2 < WqWe. IUn.ll n < C. 

II II 2 II N _ 2 

Hence, V z n — ^ Vz in L 2 (M N ). All this is enough to pass to the limit in (110.20 and 
show that the limit u* is indeed %j. □ 



10.2 <7 — y 2 

We now study the upper limit a — > 2. 
Theorem 10.2 The map S is also continuous at a = 2. 
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Proof. Without loss of generality we keep m fixed and let a n — » 2~. As before, it is 
enough to prove the convergence of the semigroup for bounded functions. The proof 
uses the extension technique introducing the vertical variable y > and considering 
the extended problem in the upper-half space ( 17.ip . The convergence, in L 1 1 oc (IR+ +1 ), 
of the sequence w n = E an (u™) to some function m;*, as well as that of Vw n to Vu>* 
in in compact sets, works as before. 

To identify the limit of the trace u* = (Tr^*)) 1 /*™, we follow an idea from [TT] . 
We take a factorized test function, £,{x)ri(y) , where t] is a cut-off function, t](y) = 1 
for y < 1, rj(y) — for y > 2, and £ is a usual test function in the x variables. 
Using that the measures {2 — a n )y l ~ Un dy are probability measures on (0, 1) converging 
(in the weak-* sense of measures) to a Dirac measure So, applied to the sequence 
¥n{y) = v(y)^xUn(x,y), we finally arrive to 

/ u£dx + [ (V<\V0cfe= / g£dx, 
Jm. N Jr n Jr n 

since linio-^2- = 1. This is the weak formulation of the equation m» — Au™ = g. □ 



11 Comments and extensions 

• Limit. The limit case a — > looks very interesting. Notice that when a = 0, the 
critical exponent is m* = 1. Formally, the limit equation for a = and m > 1 is the 
ODE 

(H.l) ^ + 1^-^ = 0. 

In case the initial datum u(x, 0) = f(x) is a function defined pointwise, the ODE can 
be explicitly solved, giving the formulas 



(11.2) u{x,t) 



(Ifix)] 1 -™ + (m - l)t)- 1 /( m - 1 )sign(/(a;)) if m > 1, 
f(x) e _t if m — 1. 



We see from these formulas that mass is not conserved; it decays instead. Thus we 
can not have convergence in L 1 (IR Ar ) as o — > + . Observe also by passing that when 
m > 1 the decay rate t -1 ^" 1-1 ) of ||it(-,£)||cc agrees with that of formula (18.61) . 

On the other hand, the question of continuous dependence below the critical expo- 
nent m*, in some weighted norm, is interesting and will be the subject of a future 
work. 

• Bounded domains. We have presented the basic facts for a theory in Sections H] 
and 17.11 We point out that there are other ways to understand the Cauchy-Dirichlet 
problem in a bounded domain with homogeneous 'boundary data'. For example, one 
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may look for solutions to the problem 




u m (y,t) - u m (x,t) 



dy, 



x e Q, t > 



x e R N \ Q, t >0, 



x G f2, 



which is different from Problem (14.1jl . This approach has been recently used by Kim 
and Lee [28], [29], who have addressed some important issues in that framework, such 
as existence, regularity and asymptotic behaviour. 

The continuous dependence proofs of the last section extend to the case of solutions 
of Problem (Oj) . 

• Related work. Recently, Cifani and Jakobsen have studied the existence of 
solutions of the diffusion-convection equation d t u + V ■ f(u) + (— A) a ^ 2 A(u) = 0, in the 
framework of (Kruzhkov-style) entropy solutions [18] . However, their assumptions on 
the nonlinearities exclude (II. ip . unless m = 1. 
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